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Abstract 



We study analytically the ground-state stability of a Bose-Einstein con- 
densate (BEC) confined in an harmonic trap with repulsive or attractive 
zero-range interaction by minimizing the energy functional of the system. In 
the case of repulsive interaction the BEC mean radius grows by increasing 
the number of bosons, instead in the case of attractive interaction the BEC 
mean radius decreases by increasing the number of bosons: to zero if the sys- 
tem is one-dimensional and to a minimum radius, with a maximum number 
of bosons, if the system is three-dimensional. 



PACS Numbers: 03.75.Fi, 05.30.Jp 
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In the last three years there has been a renewed interest in the Bose- 
Einstein condensation due to the spectacular experiments with alkali vapors 
^'^Rb, '^^Na and "^Li confined in magnetic traps and cooled down to a temper- 
ature of the order of 100 nK^^'^^'^^. Numerical studies of the Bose-Einstein 
condensate (BEC) have been performed for the ground state'^-*'^-''^-* and the 
collective excitations^-* 

In the present paper we analyze analytically the ground-state stability of 
the BEC by minimizing the energy functional with respect to the standard 
deviation of a Gaussian trial wave-function. Our analytical results are in 
good agreement with the numerical calculations of Edward and Burnett^-* 
and also of Dalfovo and Stringari^^ Wc are able to estimate the maximum 
number of bosons for which the BEC is stable. Moreover we find strong 
differences between the one-dimensional and the three-dimensional cases. 

For the alkah vapors the range of the atom-atom interaction is believed 
to be short in comparison to the typical length scale of variations of atomic 
wave functions. The atom-atom interaction can be replaced by an effective 
zero-range interaction potential^^ 

U{r - r') = B6{r - r') . (1) 

Such an effective potential leads automatically to s-wavc scattering only, 
where B — 2nh^a/m is the scattering amplitude and a is the s-wavc scat- 
tering length. This scattering length is supposed to be positive for ^'^Rb and 
'^^Na but negative for '^Li. This means that for ^"^Rb and '^^Na the inter- 
atomic interaction is repulsive while for ^Li the atom-atom interaction is 
effectively attractive^-* 
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By applying the theory of weakly interacting bosons^^-* , the Hamiltonian 
operator can be written 

H = -|^V^ + K.t(r) + Vm{v) , (2) 

where Kxi(r) is the external potential of the trap and Kn is the mean-field 
self-consistent potential, given by 

Vm{r) = j dv'\ij{v')\^U{v-v'), (3) 

where is the wave-function of the BEC. Although the actual experimen- 
tal traps are anisotropic^^'^^'^\ it is simplest to consider an isotropic harmonic 
trap (the effect of the anisotropy can be treated in perturbation theory). The 
bosons are alkali atoms in the trap 

Vext{r) = -^r^ , (4) 

with zero-range effective interaction, thus by using the equation (1) we get 

VM^B\^{v)\'. (5) 

The mean energy of the system is given by the Gross-Pitaevskii functional^^^ 

K[il:] = J dYil:*{Y)Hil:{Y) , (6) 

and we can study the ground state stability by imposing the minimization 
of the energy functional 

5K[^l;] = . (7) 

The main point of this paper is the choice of the trial wave-function 
for the energy functional. We choose a Gaussian wave-function with a free 



parameter cr, which is the standard deviation of the Gaussian, i.e. the mean 
radius of the condensate. In fact, for a — J— the test function is the 
ground-state function of the non-interacting system. Moreover, previous 
numerical results show that a Gaussian is a good approximation of the true 

ground-state function of the BEC'^-''''"^^'*'^ 

Let us start with the one-dimensional case. We choose the following trial 
wave-function 

V'(x) = Cexp(^), (8) 
with the normalization condition 

j dx \il:{x)\^ = N , (9) 

where N is the number of bosons, from which we obtain 

c" = ^- (10) 

By inserting this trial wave-function in the energy functional, after some 
simple calculations, we find 

K ^ -i )N— + -{ )N(T^ + - -7777-. (11) 

The minimum of the energy functional with respect to the standard deviation 
(7 is obtained by imposing the following condition 

^ dK ,h\^l ,mu;\^^ BN^ 1 , , 

from which we get the formula 
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It is easy to see that for any B there is only one solution of the equation 
^ — 0. The second derivative is positive when ^ = 0, namely the 
solution is stable, i.e. a minimum of the energy functional. In particular, for 
N — we have a = J— and if S > then cr — > oo for A?" — > oo, while if 
5 < then a ^ for ^ oo. 

Now we consider the three-dimensional case. We choose again a Gaussian 
trial wave-function 

V'(r) = Cexp(I^), (14) 
with the normalization condition 

^ dr IV'WI^^TV, (15) 

from which we find 



1' 



The resulting energy functional is slightly different from the one-dimensional 
one 

We find the minimum of the energy functional with respect to the mean 

radius a by imposing the following condition 

^ dK 1 ,mcj2^_ BN'^ 1 , , 

= — = -3 W— + 3 )Na-2,- — t-ttt^ , (18) 

from which we obtain the formula 

By studying the function we observe that for B > Q there is only one 
solution (intersection with the cx-axis), which is a minimum of the energy 
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functional (stable solution), instead for i? < there are two solutions: one is 
a maximum (unstable solution) and the other a minimum (stable solution) 
of the energy functional. 

It follows that for S > 0, when = we have a = J— and u — > oo for 
N ^ oo. For B < the mean radius a decreases by increasing the number 
of bosons A*" to a minimum radius cr*"*", with a maximum number iV™"^ of 
bosons. For greater values of N the condensate becomes unstable and there 
is the collapse of the wave-function. It is not difficult to obtain the critical 
number of bosons. We put B — —B and get 

from which we have the minimum radius 



5V4 

and the maximum number of bosons 



55/4 B 2m\rmj ^ ' 

Thus we obtain an analytical formula of the maximum number N"^""^ of 
bosons for which the condensate, with attractive interaction, is stable. We 
have seen that in the one-dimensional case iV™'°^ = 00 and cr""" = 0. 

For ^Li the scattering length is a = —14.5 A, the axial frequency of the 
trap is LOa/i^n) — 117 Hz and the transverse frequency is a;t/(27r) ~ 163 
Hz (see Ref. 3 for further details). We estimate the critical number of 
particle by using u;/(27r) ~ 120 Hz for the frequency of our isotropic trap: 
we find N"^"-^ ~ 1400. This result is in good agreement with the numerical 



calculations of Dalfovo and Stringari®\ It is important to observe the the 
number of particles in the BEC reported in the experimental work of Ref . 3 is 
an order of magnitude higher than our critical value. As suggested by Dalfovo 
and Stringari^\ the discrepancy between the experimental finding of Ref. 
3 and the predictions of the Gross-Pitacvskii theory could be significantly 
reduced if one assumes the existence of a vortex in the atomic cloud^^'' . 

In conclusion, we have studied in the mean-field approximation the ground- 
state stability of a gas of weakly interacting bosons in a harmonic trap with 
a zero-range interaction. To minimize the energy functional, we have cho- 
sen a Gaussian trial wave-function, where its standard deviation represents 
the mean radius of the condensate. In the case of repulsive interaction the 
mean radius of the condensate grows by increasing the number of bosons. 
In the case of attractive interaction the mean radius decreases by increasing 
the number of bosons but in different ways depending on the dimensionality. 
Our results suggest that one must be very careful when tries to simplify a 
many-body problem by reducing its dimension. 

Our technique can be applied also to bosons with non-local interaction 
and to systems with different shapes of the trapping potentiaP^^. In the 
future will be interesting to investigate analytically the collective excitation 
of the condensate, at least in the semiclassical approximation^"^)'^^^. 
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